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ON COEFFICIENTS OF POWERS OF POLYNOMIALS AND 
THEIR COMPOSITIONS OVER FINITE FIELDS 

GARY L. MULLEN, AMELA MURATOVIC-RIBIC, AND QIANG WANG 


Abstract. For any given polynomial / over the finite field Fg with degree 
at most q — 1, we associate it with a q x q matrix A{f) = (ciife) consist¬ 
ing of coefficients of its powers modulo — x for 

k = 0,1;---!?”!- This matrix has some interesting properties such as 
^ f) — ^{ 1 )^( 9 ) where {g o f)(x) = g(f(x)) is the composition of the 
polynomial g with the polynomial /. In particular, = (A(/))^ for any 

fc-th composition of / modulo x^ — x with fc > 0. As a consequence, we 
prove that the rank of A(f) gives the cardinality of the value set of /. More¬ 
over, if / is a permutation polynomial then the matrix associated with its 
inverse A(/(-i)) = A(/)-i = PA(f)P where P is an antidiagonal permuta¬ 
tion matrix. As an application, we study the period of a nonlinear congruential 
pseduorandom sequence a = {ao, ai, 02 ,...} generated by an = f^^^ao) with 
initial value ag, in terms of the order of the associated matrix. Finally we 
show that A{f) is diagonalizable in some extension field of Fq when / is a 
permutation polynomial over F^. 


1. INTRODUCTION 

Let Fg be the finite field of order q = where p is a prime number and n is 
a positive integer. Let f{x) = ^>6 a polynomial over Fg with degree at 

most q—l. To compute its composition with another polynomial g{x) = bix’', 

we can either use interpolation to obtain its expression directly, or calculate all the 
powers /(x)* (mod x® — x) in the expression {g o /)(x) = bi{f{x)y. 

Denote by 

qr-l 

(/(x))^ = ^ OifcX* mod (x”^ - x) 

i=0 

the fc-th power of the polynomial /(x) for fc = 1, 2,..., g — 1. Denote by /g the zero 
polynomial in Fq[x]. If / ^ /g we will define (/(x))° = 1 and /o(x)° = 0. 

For any polynomial /(x) = associate a coefficient vector vf with 

it, namely, 

Vf = (ao,ai,... ,0,-1)'^. 
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We define a, q x q matrix associated with f{x) ^ fo{x) by 




aoi 

002 

00,5-2 

Oo, 5 -l 

ail 

012 

01^5-2 

Ol,5-l 


05-2,2 

05-2,5-2 

05-2,5-1 

Oqr-l.l 

05-1,2 

05-1^5-2 

aq — l , q—l 


where the /c-th column consists of the coefficients of the (fc — l)-th power of f{x). In 
particular, we define ^(/o) to be the zero q x q matrix. We note that we can build 
the matrix A{f) by directly computing each of the fc-th powers of f{x) modulo 
x'^ — X. For example, finding each column of A{f) takes bit operations 

using the result of Kedlaya and Umans [5]. On the other hand, using Lagrange’s 
interpolation /(a:)^ = (l ~ can obtain the explicit 

expression for all the entries of A{f). Namely, for all 1 < i,j < g — 1, we have 

= -Ea6F, andooj = f{0y = Eaer, /(“V 

The Bell matrix of an analytic function / is an infinite matrix defined as 




j! 


dx^ 




. rc=0 


where (/(x))^ = B[f]jkX^■ It is sometimes called a Jabotinsky matrix. The 
transpose of a Bell matrix is called a Carleman matrix., which is often used in 
iteration theory to find the continuous iteration of a function [8]. 

In this paper we show that our matrix A[f ) of a polynomial / over Fg is indeed 
a finite field analogue of the Bell matrix. Some fundamental properties in terms 
of the composition of polynomials are proved similarly. Moreover, we derive a few 
results specifically related to finite field theory. In Section[2]we show that the matrix 
associated with the composition of two polynomials over a finite field is the product 
of two associated matrices. That is, A{g o /) = A{f)A{g). As a corollary, we prove 
that the value set size of any polynomial / over Fg is the rank of its associated 
matrix A{f), which is equivalent to an earlier result of Chou and Mullen [3], which 
deals with the transpose of the (I, l)-minor of A{f). In Section [3l we concentrate 
on permutation polynomials over F^. In particular, we prove that the associated 
matrix for the compositional inverse satisfies A{f^~^'>) = PA{f)P, where P 
is an antidiagonal permutation matrix defined by Pi,q-i = 1 for i = 1,... ,q and 
zero otherwise. Moreover, we show A{ f) is diagonalizable in some extension field of 
F,j. Throughout this paper, we note that /^(x) or (/(a:))^ denotes the fc-th power 
of f{x) modulo x‘^ — X, while f^^\x) denotes the fc-th composition of f{x) modulo 

x^ — X. 


2 . The matrix of a composition of polynomials 

First we derive the following obvious result. 

Proposition 1. Let f{x) = ^ 9{x) = ^ lF 5 [a;]. 

Then 

Vgof = A{f)Vg. 

Proof. The {i + I)-th coordinate of A{f)vg is given by {A{f)vg)i+i = 

On the other hand, we obtain g o f{x) = X^o ^k{f{x)y = X^o XXo = 
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J2‘i=oiJ2l=o^kaik)x\ Therefore we obtain {vgof)i+i = J2l=o^ikbk = {A{f)vg)i+i 
for * = 0,1,... g — 1. □ 

Theorem 1. Let g{x) = X]i=o G Fg[a;] and f{x) = € Fq[a;]. Let 

{g o f){x) = g(f{x)) be the composition of g with f. Then 

A{gof) = A{f)A{g). 

Proof. By Proposition [l] we see that A{f)vgk = Vgk^f for any fc-th power of the 
polynomial g. Let ak{x) = x^. Because the composition of polynomials is an 
associative operation, we have g^ o f = (ak o g) ° f = <Xk o {g ° f) = {9 ° fY- 
Therefore A{f)vgk = V(^gof)i‘ for all k = 0,1,2g — 1. Partitioning the matrix 
A{g) with columns VgO, Vg, Vg 2 ,..., Vgq-i, we derive 

A{f)A{g) = (^A{f)VgO , A{f)Vg, A{f)Vg2 A{f)Vgq-l ) 

~ ’ ''^{gof) ■> '^{g°fP ) ■ • ■ > ^(so/)<i-i ^ = A{g o /). 

□ 

We recall that /^(x) denotes the fc-th power of f{x), while denotes the 

fc-th composition of /(x). 

Corollary 1. For any given polynomial f £ F,j[x] we have that A{f^^'>) = (A(/))^, 
for any k = 1 , 2 ,.... 

This provides an algebraic way to study the composition of polynomials in terms 
of multiplication of matrices. Although the matrices associated with polynomials 
are large and costly to build, this still gives us some interesting theoretical conse¬ 
quences. We note that the transpose of the (1, l)-minor of A{f) was earlier studied 
by Chou and Mullen [3]. They gave a result on the size of the value set of / in 
terms of the rank of the (1, l)-minor of A{f); see also page 234 in [TJ]. However, 
our proof is different. 

Corollary 2. Let f be a polynomial over a finite field F^ and \Vf\ be the size of 
the value set Vf = {/(a) | a £ F^} of f. Then \Vf\ = rank{A{f)). 

Proof. If /(x) G Fq[x] is not a permutation polynomial then we define D = Vf and 
let g £ Fq[x] be a nonzero polynomial of least degree m such that g : D ^ {0}. Let 
g{x) = bmX"^ + bm-ix^~^ -I- ■ • ■ -b bix -b 6o- Then we have g o /(x) = 0, and thus 
A{f)vg = 0 by Proposition [T] This means that the first m -b 1 columns of A{f) 
are linearly dependent and thus the coefficients of g determine a linear dependence 
among the polynomials 1, /(x), /^(x),..., /"“(x) in the sense that = 

0. Moreover, (/(x))°,/(x),..., (/(x))*””^ are linearly independent because g{x) 
is the lowest degree polynomial such that g o f = 0. Therefore, rank{A{f)) = 
deg{g{x)) = \Vf\. 

11 f £¥q [x] is permutation polynomial, then all the powers of / and correspond¬ 
ing columns of A{f) are linearly independent. □ 

Corollary [5] states that the size of the value set of / is given by the rank of the 
matrix A(f). One would also like to know which elements c £¥q show up in the 
value set Vf of /, and if c shows up in the value set, how many times does it appear? 

Again, we consider the polynomial /(x) = op -b aix -b • • ■ -b aq-ix‘^~^ over F^. 
First we want to determine the number of nonzero solutions to /(x) = c. Let us 
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consider the polynomial h{x) = (ao + Oq-i — c) + aix + • • • + aq- 2 x'^~‘^. Then, by 
the Konig-Rados Theorem (Theorem 6.1 in |10jl. the number of nonzero solutions 
to /(x) = c is g — 1 — r, where r the rank of the (g — 1) x (g — 1) left circulant 
matrix 


C{h) 


Oq-l — C 

Ol 

aq-2 

«! 

02 

. ao + Oq-l - c 

02 

03 

Oi 

aq-3 

Oq-2 

(lq—4 

O-q-2 

OQ “t“ Og_l C 

aq-3 


Therefore, if c ^ /(O); then c appears in the value set Vf of / if and only if the 
rank of the matrix C{h) is less than q — 1. And the number of times that c appears 
in the value set Vf of / is r if and only if the rank of the matrix C{h) is g — 1 — r. 
If c = /(O) then c appears in the value set q — r times. 

Let k be the largest positive integer such that {1, /,..., /^} is linearly indepen¬ 
dent over Fq. Then obviously the rank rank{A{f)) > fc -|- 1. For example, let 
/ G be a polynomial of degree d, then it is obvious that 1,/,..., 
are linearly independent and thus the value set Vf has size \Vf \ > [{q — l)/d\ -I- 1. 
We note that polynomials satisfying \Vf \ = [{q — l)/(ij -I- 1 are called minimum 
value set polynomials. The classification of minimum value set polynomials is the 
subject of several papers; see [n El El lain]. Using the discussion after CorollaryEl 
we have the following. 


Corollary 3. Let f be a polynomial of degree d over the finite field Fg. Then / 
is a minimum value set polynomial if and only if rank(A(f)) = [(g — Ij/dJ -|- 1. 
That is equivalent to, {1,/,...,} is a basis which spans the space of all 
nonnegative powers of f. 

Let us consider the (1, l)-minor M{f) of A{f). If the Tth row of M{f) consists 
entirely of O’s or entirely of I’s, set h = 0. Otherwise for a nonzero *-th row of 
M{f), arrange the entries in a circle and define h to be the maximum number of 
consecutive zeros appearing in this circular arrangement. Let Lf he the maximum 
of the values of h, where the maximum is taken over all of the q — 1 rows of the 
matrix M(/). Using the linearly independence of these columns, we can derive a 
lower bound of the size of the value set Vf. 

Corollary 4. (Theorem 1, [4]j Let f be a polynomial over¥q and Lf be defined as 
above. Then \Vf\ > Lf + 2. 

In |I2) . Mullen fully classified polynomials /(x) over finite fields which commute 
with linear permutations, that is, f{bx -I- a) = 6/(x) -I- a. We note that A(bx + a) 
is an upper triangular matrix. Comparing the second column of A(bx + a)A(f) = 
A(f)A(bx + a), one can derive the following corollary. 
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Corollary 5. (Theorem 1, [T2]y) The polynomial f(x) = bo + bix + • • • + bq-ix‘^ ^ 
satisfies f{bx + a) = bf(x) + a if and only if 


( 1 ) 


bo{b-l) 

bs{l-b^-^) 


9-1 

—a + ^ ^ bto*, 

t=i 



(1 < s < 9 - 1) 


3 . Permutation polynomials 


Permutation polynomials over the field Fg under the operation of functional com¬ 
position form a group isomorphic to the symmetric group (Sq, o) with q! elements. 
There is a representation of the permutation polynomials in terms of circulant 
matrices such that its centralizer is commutative m, but here we consider the rep¬ 
resentation of / in terms of the invertible matrix A{ f). We note that the mapping 
/ —>■ A{f) is one-to-one. Hermite’s criterion (Theorem 7.4 in [10]) states that f{x) 
is permutation polynomial if and only if the coefficient Uq-i^k in the k-th power of 
f{x) is 0 for all A: = 1,2,... ,q —2 and f{x) has exactly one root in Fg, say /(e) = 0 . 
This means that all entries of the last row of A{f) are zero except Oq-i^q-i = 1. 
Indeed, f{x)'^~^ = 1 if a; ^ e and f{x)‘>~^ = 0 if x = e. Hence 


fix^-^ =Y,{l-{x-ay-^)=x'^-^ 

a^e 



(_a)^-i-V+^(l-(-a)^-i). 

a^e 


Therefore ag_i_g_i = 1. Moreover, ao,g_i = 1 if e 7 ^ 0, and ao,?-! = 0 if e = 0. 

We now consider the compositional inverse of a permutation polynomial / 

with respect to composition. Since A{g o f) = A{f)A{g) and the matrix associated 
with /^°^(x) = Id{x) = a; is the identity matrix, it is easy to see that A{f^~^'>) = 
A{f)~^. Moreover, we have 


Theorem 2. Let f be a permutation polynomial of¥q. Let P be the antidiagonal 
permutation matrix, i.e. P is defined by = 1 for i = 1,2, ...,q and zero 

otherwise. Then A{f^~^^) = {A{f))~^ = PA{f)P. 

Proof. Obviously, A{f^~^'>) = {A{f))~^. Denote the /c-th power of / and the inverse 
polynomial f^~^'> by f'^{x) = and {f^~'^Hx))'^ = respec¬ 

tively, for A: = 1, 2,..., g — 1. For any permutation polynomial g{x) = YZi=o cia;®, 
it is well known (see for example m) that its coefficients can be calculated by 
O = - EseF, s'^“^"*5(s), for i = 0, 1 ,..., q - 2 . 

For 1 < k < q — 2,hy Hermite’s criterion, the polynomial f^{x) must have degree 
at most q — 2. Therefore we have for 0 < i < g — 2 and 1 < A: < g — 2, 

= = - E = bq-^-k,q-l-^, 

seF, seF, 


i.e. 

Ojfc = for 0 < i < q - 2 and 1 < k < q - 2. 

Moreover, Oq-i^k = 0 forl<A :<9 — 2 and Oq-i^q-i = 1. In addition, aq-i^k = 
bq-i-k.o by the definition of A{f^~^^). 
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On the other hand, for any polynomial g{x) = Cis;®, it is well known that 

its coefficients can be calculated by * = 1,... ,q — 2 

and Co + Cq-i = — X^sgf, Hence we can compute 

sGF, sGF, 

for 1 < i < g - 2 and ao,,-i + Qq-i^q-i = - X^sgf, = - SsgF, = 1 - 

Because Uq-i^q-i = 1, we have ao,q-i = 0, which is equal to 60 , 9-1 by definition. 
Hence we have proven that Uik = bq-i-k,q-i-i, for all 0 < i < g — 1 and 1 < fc < 
g — 1. Since the multiplication by P on both sides reverses the order of rows and 
columns of A{f), it follows that = PA{f)P. □ 


Corollary 6 . Let f be a permutation polynomial and P be the antidiagonal per¬ 
mutation matrix as defined in Theorem\^ Then the matrix PA{f) is the inverse 
of itself. 

Proof. By Theorem [21 we have (H(/))“^ = A{f^~^'>) = PA{f)P. Therefore 
{P{A{f)f = 1. □ 


Corollary 7. Let S be a group of invertible q x q matrixes over Fg equipped with 
the operation A * B = B ■ A where B ■ A denotes the usual product of the matrices 
B and A. Denote by the permutation polynomial of degree at most q — 2 induced 
by a permutation it G Sq. Then the mapping ip : Sq ^ S given by ipiir) = H(/^) is 
a monomorphism and thus Sq is isomorphic to the subgroup A = {A(/^)|7r £ Sq} 
of the group S. 

Proof. It is easy to show that {S, = 1 =) is a group and that the mapping ip is injective. 
Now ipiiT oa) = o fa) = A{fa) ■ H(/^) = ip{Tr) * ip{a). □ 


Finally we comment on some potential applications of our results in sequence 
designs. For any permutation polynomial /, we define a nonlinear congruential 
pseudorandom sequence d = { 00 , 01 , 02 ,...} such that an = and oq is the 

initial value. The period of o is equal to the smallest k such that = 

/(*)(ao) for some i. Character sums of these sequences are studied in [SKniiniiiii 
[ 19 ]. For each initial value that is not fixed by /, we find the period of the nonlinear 
congruential pseudorandom sequence. If we take K as the least common multiple 
of all these periods, then we obtain f^^'> = id and thus A{f^^'>) = I. Conversely, 
if A(f)^ = I then the period of the nonlinear congruential sequence is a divisor of 
the order of the invertible matrix. Next we demonstrate the following two simple 
examples although they can be obtained easily without using these matrices. 

Let f{x) = cc"* be a polynomial over F, such that {m,q — 1) = 1. Then A{f) 
is a permutation array such that the only nonzero entry in column k is in (km 
(mod g — 1), fc) position where 1 < fc < g — 1. The period of a is well known, which 
is the order of m modulo g — 1. 
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Let f{x) = ax + & € Fp[a:], where a is a primitive element in Fp and b ^ 0. Then 


A{f) 


'1 b b^ 

0 a 2ab 

0 0 0 

0 0 0 


bP-2 

ip - 2)abP-^ 


0 


1 

(p - l)abP-^ 


{p - l)ap-2b 

1 


The matrix is an upper triangular matrix such that its eigenvalues are all the 
nonzero elements (a^, k = 1,... ,p — 1) in Fp and the multiplicity of 1 is 2. Hence 
the period of A is equal to p — 1, the least common multiple of orders of these 
eigenvalues. 

Computing the order of the matrix A{f) associated with a permutation polyno¬ 
mial / provides an algebraic way to find out the period of this kind of pseudorandom 
sequence, although the matrix A{f) itself is costly to build. For example, finding 
each column of A{f) takes bit operations using the result of Kedlaya and 

Umans [9]. We wonder whether we could overcome this drawback by pre-computing 
the initial matrix, or taking a sparse matrix, or diagonalizing the matrix. As an at¬ 
tempt, we end our paper with a diagonalization result of A{ f) over some extension 
field of F,. 

Theorem 3. Let f G Fq[x] be a permutation polynomial o/F^ such that the disjoint 
eycles Ci, C 2 ,..., Cfe of f have lengths Li,L 2 ,---,Lfc respectively. Let K be an 
extension field of F^ that contains all solutions of the equations x^' — 1 = 0 for 
i = 1,2,... ,k and tpi be a fixed primitive Li-th root of unity in K for each i. Then 
A(f) is diagonalizable with eigenvalues iff for i = 1,... ,k and j = 0,..., Li — 1. 

Proof. For each cycle Ci we pick a starting point (arbitrarily) and denote it by bg, so 
our cycle can be denoted by (bg, bi,..., b^.-i). For each j such that 0 < j < Li — 1, 
we can define 




if X = bt € Cp, 
0 ifx&¥q\Ci. 


Obviously, 


i.e., each pij produces an eigenvector of A{ f) with the corresponding eigenvalue 
(iff). Indeed, if x ^ Ci, then /(x) ^ Ci and so gij{x) = 0 = gij{f{x)). If 
X € Ci then x = bt for some t = 0,1,...,Li — 1. Then /(x) = bt+i (mod Li)- 
Thus gi,j{f{x)) = In this way we obtain a set 

{pij(x) -. i = 1,... ,k,j = 0,1,..., Li — 1} oi q polynomials in K[x\. For each fixed 
i, it is easy to see that {pij(x) : j = 0,..., Li — 1} is linearly independent because 
ipi is a primitive Li-th root of unity. Moreover, if * ^ z' then gij{x)gipji(x) = 0. 
Therefore the set of q polynomials {pij(x) : i = 1,... ,k,j = 0,1,..., Li — 1} is 
linearly independent. Because the size of A{f) is q and these pij(x)’s provide us 
q linearly independent eigenvectors corresponding to eigenvalue ipf, the proof is 
complete. □ 
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Remark 1. From the proof of Theorem\^ all polynomials g{x) € K[x\ such that 
g{f{x)) = Xg{x) for some A satisfy 

k Li — 1 

9{x) = XI XI 

i=l j=0 


Remark 2. Theoreml^ can be extended to non permutation polynomials such that 
either x or f{x) is in a cycle of the functional graph of f, that is, the tail length of 
any element in the functional graph is at most one. For each such a leaf d in the 
functional graph of f, we define the function 


9iAx) 


1 if X = d; 
0 if X ^ d. 


Obviously, d ^ Vf. Flence PiAfA)) = 0 = OpiAA /o’’ x € ¥g and thus 
9iAA derives an eigenvector corresponding to the eigenvalue 0. Together with the 
eigenvectors corresponding to the nodes in the cycles, we have q linearly independent 
eigenvectors and thus A{f) is diagonalizable. Flowever, in general A{f) is not 
necessarily diagonalizable in each of its extension fields. For example, let f{x) = 
+ a; + 1 € F5[j:]. Then 


A{f) 


/ 

1 

1 

1 

1 

1 

\ 


0 

1 

2 

1 

0 



0 

1 

3 

2 

0 



0 

0 

2 

2 

0 


V 

0 

0 

1 

1 

0 

/ 


It is easy to check that the rank of A(f) is 3 over F 5 . However, eigenvalues of 
A{f) oxer R are 5,1,1,0,0 and thus are 0,1,1,0,0 oxer F 5 . Hence A{f) can not 
be diagonalizable over any extension field 0 /F 5 . 


References 

[1] H. Borges and R. Conceicao, On the characterization of minimal value set polynomials, J. 
Number Theory 133 (2013), 2021-2035. 

[2] L. Carlitz, D. J. Lewis, W. H. Mills, and E. G. Straus, Polynomials over finite fields with 
minimal value sets, Mathematika 8 (1961), 121-130. 

[3] W.-S. Chou and G. L. Mullen, A note on value sets of polynomials over finite fields, preprint, 

2012. 

[4] P. Das and G. L. Mullen, Value sets of polynomials over finite fields. Finite Fields with 
Applications to Coding Theory, Cryptography and Related Areas (Oaxaca, 2001), 80-85, 
Springer, Berlin, 2002. 

[5] D. Gomez and A. Winterhof, Character sums for sequences of iterations of Dickson poly¬ 
nomials, Finite Fields and Applications, 147-151, Contemp. Math., 461, Amer. Math. Soc., 
Providence, RI, 2008. 

[6] J. Gomez-Calderson, A note on polynomials with minimal value set over finite fields, Math¬ 
ematika 35 (1988), 144-148. 

[7] J. Gomez-Calderon and D. J. Madden, Polynomials with small value set over finite fields, J. 
Number Theory 28 (1988), no. 2, 167-188. 

[8] P. Gralewicz and K. Kowalski, Continuous time evolution from iterated maps and Carleman 
linearization, Chaos Solitons Fractals 14 (2002), no. 4, 563-572. 

[9] K. S. Kedlaya and C. Umans, Fast polynomial factorization and modular composition, SIAM 
J. Comput. 40 (2011), no. 6, 1767-1802. 

[10] R. Lidl and H. Niederreiter, Finite Fields, Encyclopedia of Mathematics and Its Applications 
20, Cambridge University Press, 1997. 

[11] W. H. Mills, Polynomials with minimal value sets, Pacific J. Math 14 (1964), 225-241. 




ON COEFFICIENTS OF POWERS OF POLYNOMIALS AND THEIR COMPOSITIONS 9 


[12] G. L. Mullen, Polynomials over finite fields which commute with linear permutations^ Proc. 
Amer. Math. Soc. 84 (1982), no. 3, 315-317. 

[13] G. L. Mullen and D. Panario, Handbook of Finite Fields, CRC Press, Boca Raton, FL, 2013. 

[14] A. Muratovic-Ribic, Representation of polynomials over finite fields with circulants, Sarajevo 
J. Math. 1(13) (2005), no. 1, 21-26. 

[15] A. Muratovic-Ribic, A note on the coefficients of inverse polynomials, Finite Fields Appl. 13 
(2007) 977-980. 

[16] H. Niederreiter, Design and analysis of nonlinear pseudorandom number generators, Monte 
Garlo Simulation, A.A. Balkema Publishers, Rotterdam, 2001, 3-9. 

[17] H. Niederreiter and I. E. Shparlinski, On the distribution and lattice structure of nonlinear 
congruential pseudorandom numbers. Finite Fields Appl. 5 (1999), 246-253. 

[18] H. Niederreiter and I. E. Shparlinski, Recent advances in the theory of nonlinear pseudoran¬ 
dom number generators, Monte Carlo and Quasi-Monte Carlo Methods, 2000 (Hong Kong), 
86-102, Springer, Berlin, 2002. 

[19] H. Niederreiter and A. Winterhof, Multiplicative character sums for nonlinear recurring se¬ 
quences, Acta Arith. Ill (2004), 299-305 . 

Department of Mathematics, The Pennsylvania State University, University Park, 

PA 16802, USA 

E-mail address: inullen@math.psu.edu 

University of Sarajevo, Department of Mathematics, Zmaja od Bosne 33-35, 71000 

Sarajevo, Bosnia and Herzegovina 

E-mail address: amela@pmf.unsa.ba 

School of Mathematics and Statistics, Carleton University, 1125 Colonel By Drive, 

Ottawa, Ontario, KIS 5S6, CANADA 

E-mail address: wang@math.carleton.ca 



